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Abstract 

Complementary strands in DNA double helix show temporary fluctuational openings which are 
essential to biological functions such as transcription and replication of the genetic information. 
Such large amplitude fluctuations, known as the breathing of DNA, are generally localized and, 
microscopically, are due to the breaking of the hydrogen bonds linking the base pairs {hps). I 
apply imaginary time path integral techniques to a mesoscopic Hamiltonian which accounts for 
the helicoidal geometry of a short circular DNA molecule. The hps displacements with respect 
to the ground state are interpreted as time dependent paths whose amplitudes are consistent 
with the model potential for the hydrogen bonds. The portion of the paths configuration space 
contributing to the partition function is determined by selecting the ensemble of paths which fulfill 
the second law of thermodynamics. Computations of the thermodynamics in the denaturation 
range show the energetic advantage for the equilibrium helicoidal geometry peculiar of B-DNA. 
I discuss the interplay between twisting of the double helix and anharmonic stacking along the 
molecule backbone suggesting an interesting relation between intrinsic nonlinear character of the 
microscopic interactions and molecular topology. 
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I. HAMILTONIAN MODEL FOR DNA 



While it has been lonff recognized that the form of the DNA molecule is key to understand 
its biological function l|], a considerable amount of work has been recently devoted to 
explain the DNA dynamics in terms of the microscopic interactions at play in the double 
helix. A fully atomistic description of DNA, even of a short fragment, would represent a 
formidable computational task due to the huge number of degrees of freedom. Moreover, the 
specificity of DNA lies in its capability to store the genetic information and allow a reading 
of the latter through large amplitude motions which temporarily bring apart portions of the 
complementary strands. To describe these properties we need mesoscopic models at the scale 
of the base pair, the fundamental entity in the nucleotide which encodes the information. 
The one-dimensional Dauxois-Peyrard-Bishop (DPB) model has provided a fundamental 
tool for the biophysicists working in the field. 

The DPB Hamiltonian for a system of bps, with reduced mass /i, assumes the pair 
mates separation ?/„ (for the n-th base pair) with respect to the ground state position as the 
relevant degree of freedom. The inter-base pair interactions are modeled by a Morse potential 
Vjud/n) whereas the intra-base pair stacking along the molecule backbone is described by an 
anharmonic potential Vs{yn,yn-i)- 

Recently, I have proposed to apply the path integral method to a modified DPB 
Hamiltonian which includes a twist angle 6 between adjacent bases, n and n — 1, along 
the DNA backbone [j] as shown in Fig. [H Twisting is described by the angle that the bps 
rotate around the molecule axis. B-DNA at room temperature has a helix repeat of ~ 35 A 
hosting h 10 bps, hence the equilibrium twist angle is d^q = 2Tx/h ~ O.Qrad. Taking a 
short fragment, = 100, the (integer) equilibrium twist is {Tw)eq = N/h = 10 The 
Hamiltonian reads: 



TV r .9 



n=l 



+ Vsiyn, yn-l) + VM{yn) + Vsoliyn 



K 
~2 



1 + p exp [-«(?/„ + yn-i)] {yl - 2ynyn-i COS + yl_^) 
Vuiyn) = Dn{exp{-anyn) - l)^ 

VsoiiVn) = -Dnfs{tiiiah{yn/ls) - l) • (1) 
Dn and a„ are the pair dissociation energy and the inverse length setting the hydrogen 




FIG. 1: (Color online) (a) Fixed planes picture for the right-handed helicoidal model. The blue 
filled circles denote the pointlike base pairs stacked along the molecule axis with twist 0. The 
radial coordinate y„. describes the n — th base pair displacement from the ground state. The 
dashed vertical axis corresponds to the Un = configuration, the minimum for the one-coordinate 
potential VMiUn) + ysoiiUn) in Eq. ([T|). (b) The helix plane seen from above. 

bond potential range for the n-th base pair. K is the harmonic stacking whereas p and a are 
the anharmonic stacking parameters which are taken independent of the type of base at the 
n and n—1 sites. The homogeneity assumption for the stacking relies on the observation that 
both types of bps contain a purine plus a pyrimidine, the former being larger and heavier. 
Thus the AT- and GC- bps are comparable in size and weight. For the choice of the model 

n 

potential parameters see Refs.[J]. Eq. ([T]) also introduces the solvent potential Vsoi which 
adds to Vm thus enhancing by /s-D„ the height of the energy barrier above which the base 
pair dissociates. The factor = 0.3 mimics the effect of a high salt concentration which 
screens the negatively charged phosphate groups. The length 1^ = 3A defines the range 
beyond which the Morse plateau is recovered and Z)„ returns to be the fundamental energy 
scale. For ?/„ > Ig, the two strands are apart from each other and the hydrogen bond with 
the solvent is established. 
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II. PATH INTEGRAL METHOD 

The imaginary time path integral method jsl is apphed to Eq. ([1]) by introducing the 
idea that the Un can be described by paths x(rj), the latter being periodic functions of the 
imaginary time Tj, x(rj) = x(rj + /3) with /3 being the inverse temperature,. The index 
i numbers the hps along the r-axis. In fact, there are + 1 base pairs in Eq. ([1]) but 
the presence of an extra base pair uq is remedied by taking periodic boundary conditions. 
Ho = Vn, which close the finite chain into a loop. This condition is incorporated in the 
path integral description as the path is a closed trajectory, a;(0) = x{(3). Hence a molecule 
configuration is given by paths and, in the discrete time lattice, the separation between 
nearest neighbors bps is Ar = P/N. Then, Eq. ([T]) transforms onto the time axis by 
mapping: ?/„ — x{Ti) and — )■ x(rj — Ar). The paths can be expanded in Fourier series 
with cutoff Mp 



Mf 

x{Ti) = xo + ^ cos(2m7rri/ f3) + 6^ sin(2m7rri/ (3) , (2) 



m=l 



and this introduces the following physical picture: 

a) given a set of coefficients {xq, a^, bm}, the bps are represented by the configuration 
{xin),t= 1,.., N}. 

b) A set of coefficients corresponds to a point in the path configuration space thus, 
sampling the latter amounts to build an ensemble of distinct configurations for the system. 
As this is done for any temperature, we have a tool to describe the base pair thermal 
fluctuations around the equilibrium (x(rj) ~ 0). 

c) In principle the configurations ensemble for the DNA fragment is infinite as it may in- 
clude any possible combination of Fourier coefficients. For practical purposes some physical 
criteria intervene to select computationally the path coefficients defining a molecule config- 
uration and contributing to the partition function. This poses a restriction on the ensemble 
size. 

Such criteria are of two types: First, Vm excludes too negative base pair stretchings due 
to the hard core which mimics the repulsion between negatively charged sugar-phosphate 
groups. Second, the ensemble of paths has to be consistent with the thermodynamics laws. 
This means that the numerical code selects, at any temperature, a path ensemble and 
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evaluates the entropy of the DNA fragment. If the entropy is growing versus T, the code 
proceeds to the next temperature step otherwise a new partition is performed in the Fourier 
coefficients integration, a new path ensemble is selected and the entropy is recalculated. This 
is done at any T until the macroscopic constraint of the second law of thermodynamics is 
fulfilled throughout the whole investigated temperature range. I emphasize that the method 
does not put any constraint on the shape of the entropy versus T- plot aside from the 
requirement that the entropy derivative has to be positive. It follows that the path ensemble 
is a dynamical object accounting for the manyfold of molecule configurations which enter 
the thermodynamical calculation. The size of the ensemble is a measure of the cooperativity 
degree of the system. By increasing T, some bps may open and cooperatively lead to bubble 
formation along segments of the double helix. Accordingly the ensemble size is expected to 
grow versus T. 

Applying the mapping technique to Eq. ([1]), the classical partition function for the DNA 
molecule in the solvent is 



j Dxexp -f3j2[^i^{nf + Vs{x{Ti),x{ri - At)) + Vjv/(x(ri)) + Vsoi{x{t,)) 



Mp 



Dx = J dxo (-p-) j dam j dbm , (3) 



where = ^JlTf^K is the thermal wavelength. From Eq. ([3]), I compute the ensemble 
average for the displacement of the i — th base pair, < x{Ti) >, which permits to get the 
fraction / of open bps. As the UV signal changes quite abruptly when the bps dissociate, / 
is defined in terms of the Heaviside function as: / = A^~^ J2iLi'^{^ x{Ti) > —()■ 

The threshold ( yields a criterion to establish whether an average base pair displacement 
is open, < x{Ti) >> C, or not. 



III. ANHARMONIC STACKING VERSUS TWIST 

The computational method is here applied to investigate the interplay between twisting 
and stacking anharmonicity also in view of the special role given to the latter in the DPB 
model where a finite p induces a sharp denaturation transition driven by sizeable melting 
entropy [2,]. In Fig. [2]the melting profiles are plotted both for the DPB model with zero twist 
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FIG. 2: (Color online) Fractions of average displacements larger than C = 0.6 (circles), 0.8 
(squares), 1 A (diamonds) versus temperature in the DPB ladder model ((a) - (c)) and in the 
equilibrium twist conformation ((d) - (f)). Three anharmonic stacking p are assumed both in 
the left and in the right panel. In (e) and (f), p is smaller by a factor 100 than in (b) and (c) 
respectively. 

(left panel) and for the equilibrium twist conformation (right panel). Several p values are 
assumed in both cases. In the Tw = system, even very large p produce scant variations in 
the denaturation patterns pointing to a substantial irrelevance of the p driven anharmonicity 
for the DPB ladder model. Quite different is the physical picture emerging from the right 
panel regarding the {Tw)eq conformation: even slight enhancements over the p = 2 value i) 
shift upwards along T the opening of the average base pair displacements and ii) flatten the 
melting profiles suggesting that the denaturation becomes more gradual. Note that p values 
in Figs. [2t^e),(f) are two orders of magnitude smaller than in Figs. [2]|^b),(c) respectively. The 
p = 8 plot says that, even at T ~ 4007^, there are no average displacements larger than lA 
while only 35% are larger than O.gA. Thus it is the twisting that switches on the p effect. 
The latter induces those cooperative interactions along the molecule backbone which are 
peculiar of the fiuctuational openings. Accordingly, the anharmonic stacking renders the 
double helix flexible hence it does increase the molecule resilience against the whole thermal 
disruption of the hydrogen bonds. In this sense, anharmonicity is a stabilizing factor for 
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the double helix. Modehng DNA by path integrals offers the advantage to include in the 
computation a great number of molecule configurations meanwhile accounting for those 
fiuctuational effects which are key to the molecule dynamics, mainly in short fragments. 
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